NASA/TP— 2000-2 1 0387 



Thermodynamic Cycle Analysis of 
Magnetohydrodynamic— Bypass 
Hypersonic Airbreathing Engines 

R.J. Litchford and J. W. Cole 

Marshall Space Flight Center, Marshall Space Flight Center, Alabama 
V.A. Bityurin 

Institute of High Temperatures (IVTAN), Russian Academy of Science, Moscow, Russia 
J. T. Lineberry 

LyTEC, LLC, Tullahoma, Tennessee 


National Aeronautics and 
Space Administration 

Marshall Space Flight Center • MSFC, Alabama 35812 


July 2000 


Acknowledgments 


The authors of this work acknowledge the support of the National Aeronautics and Space Administration (NASA), Space 
Transportation Directorate (STD), George C. Marshall Space Flight Center (MSFC). The results of ‘h* ^arch >«jP«tod 
herein were obtained at the request of the Space Transportation Directorate through the corroborative efforts of individuals 
associated with the Advanced Space Transportation Program Office (TD15/ASTP), MSFC; the Propulsion Research Center 
(TD40/PRC), MSFC; the Institute of High Temperatures (IVTAN), Russian Academy of Sciences > Moscow, Russra, and, 
LyTEC, LLC., Tullahoma, TN. The NASA Principal Investigator was Ron J. Litchford, TD40/PKC. 

This technical publication has been reviewed and is approved. 


George R. Schmidt 
Deputy Manager 
TD40/PRC 

Marshall Space Flight Center 


Gary M. Lyles 

Manager 

TD15/ASTP 

Marshall Space Flight Center 


Available from: 


NASA Center for AeroSpace Information 
7121 Standard Drive 
Hanover, MD 21076-1320 
(301) 621-0390 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22161 
(703) 487-4650 


11 



TABLE OF CONTENTS 


1. INTRODUCTION j 

2. THERMODYNAMIC CYCLE ANALYSIS 3 

2. 1 Thermodynamic Cycle Description 3 

2.2 Thermal Design Constraints 4 

2.3 Basic Definitions 3 

2.4 First Law Considerations 6 

2.5 Second Law Considerations 9 

2.6 System Thrust Characteristics 17 

2.7 Ramjet Mode Mach Number Constraint 19 

3. REPRESENTATIVE PERFORMANCE CALCULATIONS 20 

3.1 Specific Thrust 20 

3.2 Ramjet Mode Envelope 22 

3.3 Stream Thrust Analysis 23 

4. MAJOR TECHNICAL ISSUES 24 

4. 1 Air Ionization Considerations 24 

4.2 MHD Technology Considerations 28 

5. CONCLUSIONS 30 

REFERENCES 


iii 



LIST OF FIGURES 


Generic configuration of an airbreathing hypersonic engine with MHD-energy 
bypass system. Reference stations and process terminology are indicated. Also 
shown is the enthalpy-entropy process diagram for the generic MHD-bypass 
hypersonic airbreathing engine configuration 

Specific thrust characteristics as a function of flight Mach number for 
rj N{g) = 0, 0.25, and 0.5 assuming l s (g) =r ls(a) =Q ^ and ^O.lin^ 3000 K 

Maximum flight Mach number for maintaining subsonic combustion conditions 
as a function of the generator enthalpy extraction ratio. Assumes a maximum 
static burner entry temperature of T’ 3 i im = 1 600 K 

Naturally occurring interaction parameter behind a normal shock as a function of 
flight Mach number and altitude. Based on 1976 standard atmosphere and 
equilibrium air with molecular dissociation 

Interaction parameter behind a normal shock as a function of flight Mach number 
for both unseeded and cesium-seeded equilibrium air. Calculations were performed 
for fixed altitudes of 100,000 and 125,000 ft 

Hall parameter behind a normal shock as a function of flight Mach number for 
equilibrium air. Calculations were performed for fixed altitudes of 100,000 and 
125,000 ft with B- \ 



LIST OF ACRONYMS 


emf 

electromotive force 

HTSC 

high-temperature superconductors 

MHD 

magnetohydrodynamic 

TSFC 

thrust-specific fuel consumption 


0 

P 

Y 

A 

e 

% 

% 

% 

K 

n 

p 

a 

x 

X 

¥ 

A 

a 

B 

b 


ABBREVIATIONS AND SYMBOLS 

stagnation condition, subscript 
Hall parameter 
specific heat ratio 
relative change 
electromotive force 
combustion efficiency 
enthalpy extraction/addition ratio 
isentropic efficiency 

exit-to-entrance stagnation pressure ratio 
global stagnation pressure ratio 
density 

electrical conductivity 
exit-to-entrance stagnation temperature ratio 
fraction of generator power diverted to pre-ionizer 
exit-to-entrance static temperature ratio 
cross-sectional area 

atmospheric condition; air; MHD accelerator; constant, subscript 

magnetic field intensity 

burner, subscript 

constant pressure specific heat 


v 



ABBREVIATIONS AND SYMBOLS (Continued) 


d 

diffuser, subscript 

e 

expansion, subscript 

ent 

entrance condition, subscript 

F 

thrust 

f 

fuel-to-air mass flow ratio 

f 

fuel, subscript 

8 

MHD generator, subscript 

h 

enthalpy 

1 sp 

fuel specific impulse 

j 

current density 

L 

length 

lim 

limiting condition, subscript 

m 

mass flow rate 

M 

Mach number 

n 

nozzle, subscript 

P 

pressure 

P 

pre-ionizer, subscript 

Q 

MHD interaction parameter 

If 

fuel heating value 

R 

gas constant 

Sa 

stream thrust function 

T 

temperature 

u 

velocity 

vv 

weight flow rate 

X 

reference condition, subscript 


VI 



TECHNICAL PUBLICATION 


THERMODYNAMIC CYCLE ANALYSIS OF MAGNETOHYDRODYNAMIC-BYPASS 
HYPERSONIC AIRBREATHING ENGINES 

1. INTRODUCTION 


Established analyses of conventional ramjet/scramjet performance characteristics indicate that a 
considerable decrease in efficiency can be expected at off-design flight conditions. This can be explained, 
in large part, by the deterioration of intake mass flow and limited inlet compression at low flight speeds and 
by the onset of thrust degradation effects associated with increased burner entry temperature at high flight 
speeds. In combination, these effects tend to impose lower and upper Mach number limits for practical 
flight. It has been noted, however, that magnetohydrodynamic (MHD) energy-management techniques 
represent a possible means for extending the flight Mach number envelope of conventional engines.' - ^ By 
transferring enthalpy between different stages of the engine cycle, it appears that the onset of thrust 
degradation may be delayed to higher flight speeds. Obviously, the introduction of additional process 
inefficiencies is inevitable with this approach, but it is believed that these losses are more than compensated 
through optimization of the combustion process. 

The fundamental idea is to use MHD energy conversion processes to extract and bypass a portion 
of the intake kinetic energy around the burner. We refer to this general class of propulsion system as an 
MHD-bypass engine. In its generic configuration, an MHD generator is placed between the inlet and the 
burner as a means of converting flow kinetic energy into electrical power. In this way, the overall static 
temperature rise associated with inlet flow deceleration can be actively constrained, and the propulsion 
system is able to reach a higher freestream Mach number (i.e., higher inlet stagnation enthalpy) before the 
burner entry temperature exceeds the design limit. Furthermore, given a fixed limit for burner entry 
temperature, the inlet MHD generator can decelerate the flow to a lower velocity than that attainable using 
a simple adiabatic compression process. Thus, it is possible to increase the freestream Mach number for 
which the flow remains subsonic throughout the burner, and the ramjet mode can be sustained to much 
higher flight Mach numbers. 


In this scheme, the inlet air must be ionized and have high electrical conductivity to attain effective 
MHD interaction in the generator. At moderate hypersonic Mach numbers of 12-18, it is possible to obtain 
the required levels of conductivity through equilibrium ionization of seeded air. However, at low hyper- 
sonic Mach numbers of 5-12, a fraction of the bypassed power must be diverted for operation of a 
nonequilibrium pre-ionizer. Therefore, inlet air ionization is a key technical issue with respect to practical 
viability. 



Although the MHD-bypass engine concept is of considerable contemporary interest to the 
hypersonics community, a review of the open literature has revealed only a few analytical investigations 
aimed at performance prediction of these systems . 4 " 9 These analyses demonstrate favorable performance 
characteristics under highly specific flight conditions, but a simplified thermodynamic cycle analysis wou 
also be useful for the purpose of assessing systems performance potential on a more generalized basis. 


This technical publication provides a brief description of a recently reported analysis assessing the 
performance potential and scientific feasibility of MHD-bypass hypersonic airbreathing engines using 
fundamental thermodynamic principles.' 0 The cycle analysis, based on a thermally and caloncally perfect 
gas, incorporates strategically placed magnetohydrodynamic devices in the flow path and accounts for 
aerodynamic losses and thermodynamic process efficiencies in the various engine components, 
device performance is completely described in terms of an enthalpy extraction/addition parameter and an 
isentropic efficiency. A provision is also made for diverting a fraction of the bypassed electrical power to an 
air pre-ionizer located at the inlet of the diffuser. The power consumed by the pre-ionizer is considered 
only as a variable parameter. The detailed ionization kinetics are not addressed in this paper. The analy sis, 
while certainly not applicable to detailed performance analysis, is fundamental and is suitable as an indica- 
tor of potential performance gains associated with MHD energy bypass. It also reveals the flight ac 
number range over which the system can effectively operate and suggests the range of component 
efficiencies needed for successful implementation. 


The prime technical objective of this work was to find simple closed-form solutions for the 
performance of MHD-bypass airbreathing engines and use them to expose important trends and sensitivi- 
ties as well as to establish thermodynamic feasibility, without recourse to elaborate thermochemical 
calculations. From this standpoint, our methodology adheres to the analysis philosophy and spirit of Builder s 

pioneering approach . 1 1 
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2. THERMODYNAMIC CYCLE ANALYSIS 


2.1 Thermodynamic Cycle Description 

We consider the thermodynamic cycle of a hypersonic airbreathing engine of the ramjet/scramjet 
class that is augmented with an MHD energy management system. The geometrical configuration investi- 
gated here follows the generic MHD-bypass engine concept. The various engine components are shown in 
fig. 1 along with selected reference stations at critical axial positions along the engine flowpath. This figure 
also includes the entropy-enthalpy diagram for a modified Brayton cycle as represented by a sequence of 
eight process trajectories. 

Hypersonic Airbreathing Engine with MHD Energy Bypass System 


Power 

Conditioning 



12 3 4 5 6 10 

MHD MHD 

Diffuser Generator Combustor Accelerator Nozzle 



Figure 1 . Generic configuration of a hypersonic airbreathing engine with MHD-energy bypass 
system. Reference stations and process terminology are indicated. Also shown is the 
enthalpy-entropy cycle diagram for the generic MHD-bypass hypersonic airbreathing 
engine configuration. 
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These process trajectories may be summarized as follows: (a->l) air pre-ionization and heat 
addition (1— >2) adiabatic compression and flow deceleration, (2 >3) MHD conversion of total enthalpy to 
electrical power and deceleration of the flow, (3->4) constant static pressure and frictionless heat addition, 
( 4 _» 5 ) adiabatic expansion to prevent the temperature from exceeding design limits in the accelerator, 
( 5 _» 6 ) MHD conversion of electrical power to total flow enthalpy and acceleration of the flow, and (6— >10) 
adiabatic expansion and acceleration of the exhaust flow. 


2.2 Thermal Design Constraints 

Presently, materials used for the walls of combustion chambers and nozzles cannot tolerate 
temperatures much above 1200 K. However, unlike gas turbine engines, the surfaces of ramjet/scramjet 
combustors can be kept much cooler than the main fluid stream by providing a shielding layer of relatively 
cool air near the wall. Therefore, this class of engine can accept higher peak burner temperatures and can 
operate at higher flight Mach numbers. In this way, the relative performance and operating range of ramjet/ 
scramjet engines is greatly improved and extended over the gas turbine engine. As the flight Mach number 
continues to increase, however, the combustor inlet temperature also increases until, at some limiting Mach 
number, the peak cycle temperature begins to approach the temperature limit set by the wall materials and 
cooling methods. Because the heat energy added by fuel combustion can generate unacceptable temperature 
levels at the burner exit, we introduce a temperature design constraint on the peak combustor 
stagnation temperature in the form: 


To, 4 - ^O.lim {l) 

Because the stagnation temperature is always increased in the MHD accelerator and this device is 
subject to similar thermal design limits as the burner, it is possible to introduce an additional design 
constraint on the peak accelerator temperature in the form: 

%^ 0( iim (2) 

Even if material temperature limits could be extended, the static temperature at the burner entrance 
cannot be increased indefinitely. At temperatures approaching 2000 K, energy losses due to unequilibrated 
dissociation begin to impair performance and eventually overwhelm any benefits associated with increased 
cycle temperature. Based on the well-known thermodynamic characteristics of air, the maximum allow- 
able burner entry temperature normally ranges from 1400-1700 K. 12 Therefore, a typical design value is 
= 1600 K. For analysis purposes, we impose this practical thermodynamic constraint on design by requiring 
that the static temperature at the burner entrance not exceed a specified limiting value: 


h ^ ^lim ^ 

This limiting value should be low enough that air may be approximated as a thermally perfect gas 
with no dissociation effects during the entire inlet conditioning process. 


4 



2.3 Basic Definitions 


To simplify the analysis, we assume that the working medium can be treated as a pure substance 
throughout the engine flow path and that the working medium is thermally and calorically perfect. Because 
practical aerodynamic (non-isentropic) losses can lead to significant entropy increases and stagnation pres- 
sure losses along the flow path, it is necessary to introduce various process efficiencies into the analysis. 
These efficiencies are introduced in terms of a stagnation pressure ratio (n= p 0 exil /p 0 in ) for each engine 
process. Under these assumptions, it is possible to develop relations describing the thermodynamic process 
trajectories in each component of the engine and to reach a simple closed-form solution for fuel specific 
impulse (I sp ) and thrust-specific fuel consumption (TSFC). 

We begin the development by specifying the freestream total temperature and total pressure in the 
stream tube of the engine inlet. These properties can be defined in terms of static conditions, which are 
altitude dependent, and the flight Mach number M a using the following well-known gasdynamic 
relationships: 


T 0,a = T a 


1 + — — -M, 


*) 


(4) 


and 


PO,a ~ Pa 




(5) 


where y is the specific heat ratio. The objective is to track the variation in stagnation properties through the 
various engine components using the first and second laws of thermodynamics. 

Before proceeding with the analysis, it is useful to note that any MHD device can be described in 
terms of two fundamental thermodynamic parameters: an enthalpy extraction/addition ratio (%) and an 
isentropic efficiency {r\ s ) which accounts for aerodynamic losses as well as joule-dissipation effects. The 
enthalpy extraction/addition ratio is defined as the ratio of the stagnation enthalpy change in the device to 
the entrance stagnation enthalpy: 


VN -7 

”0,ent 


( 6 ) 


The isentropic efficiency represents the degree to which the actual process approaches an 
isentropic process. For a generator (energy extraction), it is defined as the ratio of the actual change in 
stagnation enthalpy to the ideal (isentropic) change in stagnation enthalpy that would accompany the same 
change in pressure: 


^s(g) ~ 


^khactua[ 

Mo 

, isentropic 


(7) 
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For an accelerator (energy insertion), it is defined as the ratio of the ideal to the actual. 


Vs(a) 


^OJsentropic 

actual 


( 8 ) 


Note that the enthalpy extraction ratio and isentropic efficiency performance parameters fully 
define the operational characteristics of the generator and accelerator, including joule-dissipation losses. 
Representative values can be specified based on historical data from past MHD research. 


2.4 First Law Considerations 

The pre-ionizer takes a fraction (%) of the energy produced by the MHD generator and introduces it 
to the inlet flow through an external energy addition process. The purpose is to increase the ionization level 
and electrical conductivity of the air before it enters the MHD generator. In the present formulation, we 
assume that no compression occurs in the pre-ionization region. Application of the first law of thermodynamics 
to the inlet stream tube in the pre-ionization region yields an energy balance in the form. 


= Mafaa + XVNlg)^) ^ 

where h 0 is the specific stagnation enthalpy, m a is the mass flow rate of air through the engine, and % (g) 
is the enthalpy extraction ratio for the MHD generator. The quantity %( g) *o,2 is the s P eclfic electrical 
energy extracted by the MHD generator. Since the flow is adiabatic in the diffuser, h 02 =h 0] , we can obtain 

the form: 


Tqj _ fr),i __ 1 

T 0,a fy),a 1 ~ XV. N(g ) 


( 10 ) 


The compression process in the diffuser is assumed to be adiabatic such that the stagnation 
temperature remains invariant: 


Tq,2 ~ T 0,\ ~ 


l 0,a 


XV N (g) 


(ID 


By extracting flow enthalpy and converting it to electrical power, the MHD generator acts as an 
inlet flow conditioner for the engine. In the MHD generator, both total enthalpy and total pressure are 
decreased. Application of the first law of thermodynamics to the MHD generator yields an energy balance 

in the form: 


"i a fy),3 = rh a[ h 0a~VN(g) h O.2) 


( 12 ) 
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where %( g) ^ 0 ,2 1S the total enthalpy converted to electrical power per unit mass flow of air through the 
generator. In terms of stagnation temperatures, eq. (12) can be written in the form: 


r 0,3 - 70 , 2 f 1 - tyy(*)) = T 0 ,a 

\ \g)j i -xn N ( g ) (I3 > 

Where we have introduced eq. (11) to obtain a result in terms of the freestream stagnation 
temperature. 

In the burner, the stagnation temperature is increased as the chemical energy of the fuel is released 
as combustion heat. This process is assumed to occur at constant static pressure according to the assumed 
thermodynamic cycle although actual implementation might correspond more closely to burning in a 
constant-area duct. The energy equation applied to this idealized combustion process, neglecting the 
enthalpy of the incoming fuel, is 


(m a + m f )ho A = rh a h $2 + r} c ihfq f 


( 1 + /> 0,4 =* 0 , 3 +*? cfaf 


0 + f) C p T 0 A = CpT[ 0,3 + r\ c fy f 


(14) 


where mj represents the mass flow rate of fuel, qj is the fuel heating value, Tj c is the combustion efficiency 
of the burner, C p is the constant pressure specific heat of the working fluid, and/is the fuel-to-air mass flow 
ratio. The combustion efficiency is included to account for the fact that some of the chemical energy of the 
fuel may not be released due to inadequate mixing or reaction time. The parameter may also be used to 
reflect energy losses associated with chemical dissociation. For constant specific heat, eq. (14) can be 
solved for the stagnation temperature at the burner exit: 


T _ r <u +M//C,) _ 

'0,4 — r 0,li 


(15) 


Alternatively, eq. (14) can be solved for/in the form: 


/ = 


(Toa/Tq, 3 )-\ 

(Vc<!f/C p T 0t3 )-T 0A /T 0 ' 3 


( 16 ) 
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where T n , can be obtained from eq. (13). It should be apparent that the inequality of eq. (15) can be 
satisfied by limiting the value of T 0>3 through the extraction of enthalpy in the MHD generator and/or 
constraining the fuel-to-air ratio /such that 




1 (^0,3/^O.lim) 

{jlcQ / /CpTo, lim ) — 1 


( 17 ) 


A weak expansion component may be included after the burner to expand and cool the flow enough 
such that, after passage through the MHD accelerator, the static temperature remains below specified 
design limits. We assume that the expansion is adiabatic such that the stagnation temperature remains 

constant (Tqs=Tq^). 

The MHD accelerator provides a mechanism for enthalpy addition in which the kinetic energy 
extracted from the inlet is recovered as an augmentation to overall engine thrust. In the accelerator, both 
total enthalpy and total pressure are increased. 

Application of the first law of thermodynamics to the MHD accelerator yields an energy balance in 
the form: 

(rh a + rhf)hQ 6 = ( m a + m/)[/io,5 + ^(a)^] 

7b,6=%( 1 ' +%(<.)) (18) 


where we have introduced the enthalpy addition ratio (% (a) ). The product ri N{a) h 05 is the total enthalpy 
added to the flow per unit mass flow rate of combustion gases through the accelerator. 

Because the electrical energy available for acceleration equals the total power produced by the 
generator minus that diverted for pre-ionization, % (fl) is directly related to the enthalpy extraction ratio of 
the generator according to the following power balance: 

n N (a)(< + "'f) h 0,5 = 0 - X)n N{g )rnahoa ( 19 > 


or 


%(a) r 0,5 


j^ l -x)n N {gf 0,1 


( 20 ) 
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to^taato 6 haVe ** f3C ‘ lh “ 7 ’°’ 2=7 ’ 0J ' UP °" Subs,i "" io " of 1 «>■ < 20 > im ° eq. (18) we eliminate 


T 0,6 ~ Tq,5 + 


' 1 > 


vi + /y 




Eliminating 7 0>1 using eq. (10) and noting that T 05 =T 04 , we obtain the relation 


( 21 ) 


2{),6 - 2 q ,4 + 


1 f - xK {g) 

1 + /J 1 - OTjvU ) 


^0,a 


( 22 ) 


The expansion process in the nozzle is assumed to be adiabatic such that the stagnation temperature 
remains invariant (7^ 10 =7’ 0 6 ). 


2.5 Second Law Considerations 

As the propellant flows through the engine, irreversibilities result in an increase in entropy and a 
reduction in stagnation pressure. Therefore, real aerodynamic losses (or gains) may be accounted for through 
the introduction of a stagnation pressure ratio (n= p 0jait /p 0M ) for each engine process. Expressions for 
t ese n parameters can be developed from fundamental thermodynamic and gasdynamic principles as 
demonstrated for conventional engine components by Heiser and Pratt 12 and for MHD devices by Bityurin 
et al. Because the n parameters for the MHD devices in the engine flow path introduce some novel 
considerations; however, we briefly develop the appropriate expressions for these components below. 

2.5.1 Diffuser 


The compression or diffuser efficiency rf s{d) is defined as 


tfs(d) ~ ' 


(AA), 


isentropic _ ^2 ~ h x 
actual ^2 - 


<1 


(23) 


where h x is the enthalpy at the inlet necessary to achieve the same static pressure rise in an isentropic 
compression process. Equation (23) may also be written in the form: 


Vs(d) = 




Yd~ 1 


(24) 
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or 


^ = Vd(l~ r ls(d)) + r ls(d) 

l \ 


( 25 ) 


Here, we have introduced the diffuser static temperature ratio \ff d =T 2 /T X as a new independent 
variable. 

Applying Gibbs equation for the idealized isentropic compression process yields 


p R T U-lJ T 


Thus, integration of eq. (26) through the diffuser gives 


P 2 _ _ 
P 1 


f 'T' \ 

h_ 

\Jx J 


r 

y - i 



y 

y - 1 


(26) 


(27) 


Equations (25) and (27) can now be combined to eliminate the ratio T X IT X . The result is an 
expression for the diffuser static pressure ratio. 


Pl_ — 
P\ 


Vd 


Vd^-Vs^ + TlsV) 


y-1 


(28) 


Now, by using the Mach number relation between the stagnation and static pressures, we can write 
the stagnation pressure ratio for the diffuser in the form. 


*d 


P 0,2 _ P 2 
Po,\ P\ 


1 + 


1 + 


Izl 

2 

y-1 


2 



y 

y-1 



(29) 
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Since the total temperature is conserved, we can also use the Mach number relation between the 
stagnation and static temperatures to obtain the following expression for y/^: 



1 + 


1 + 


r-i 

2 

7~1 

2 



Combining eq. (29) and eq. (30) gives 


(30) 




y 

_ Pi ( 1 V- 1 


(31) 


Finally, eq. (28) can be substituted into eq. (31) to obtain the desired result for the stagnation 
pressure ratio in the diffuser: 


+ (32) 
where we now indicate a limiting value for the static temperature ratio. 

2.5.2 MHD Generator 

Con sider the MHD generator in which the stagnation pressure must decrease. We begin with the 
Gibbs equation for an ideal (isentropic) process as applied to stagnation properties 

dpp Y a'/ip 

PO 7-1 (33) 

and integrate through the generator to obtain the stagnation pressure ratio in terms of the total enthalpy 
ratio: 




P 0,3 

P0,2 


fy),3 

^ 0,2 


7-1 


(34) 
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Utilizing the definition of rf N{g) given by eq. (6) and the definition of 7J s(g) given by eq. (7), we can 
form the following relationship between generator performance parameters: 

^lN(g) (^0,2 *0,s\ isentropic (^0.2 *0^) actual 
Tl s(g) (*0,2 '^actual *°’ 2 

( ^ 0,2 ~ ^ 0,3 
v *0,2 Jisentropic 


= 1 - 


^ fy),3 ^ 

\* 0,2 ) 


(35) 


isen/ro/nc 


Upon substituting eq. (35) into eq. (34), we obtain the desired expression for the stagnation 
pressure ratio in the MHD generator: 


*g = 


P0 , 3 _ 

P 0,2 


j _ V -1 

^s(g) J 


(36) 


2.5.3 Burner 

Because of the heat addition occurring in the burner, there is an increase in entropy and a 
corresponding loss in total pressure. This drop in total pressure is a direct manifestation of ay eig 
heating loss which obeys the following differential relationship: 


4po = _Xa/ 2 ^- 

Po 2 To 


(37) 


For constant pressure heating in the burner, the Rayleigh heating loss law leads to the following 
expression for the stagnation pressure ratio in the burner 12 (pp. 77-80 of the reference). 


P 0,4 

n b = = 

P 0,3 


i + ^V 


1-— 

n. 


-y 

r-i 


(38) 


where r h =T 0 JT 0 3 is the stagnation temperature ratio in the burner. The value of the burner entry Mach 
number Af, depends upon the design configuration of the MHD generator. 
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2.5.4 Cooling Expansion 


A weak expansion component may be included after the burner to expand and cool the flow enough 
such that, after passage through the MHD accelerator, the static temperature remains below specified 

^fve'™do^>T„!ty ® ree ° f eXPa " Si0n W ° Uld VCry ' OW - WC aMic ‘P a,e that *. wuld la ^ e on a 

Assuming adiabatic expansion, the stagnation temperature remains constant throughout the nozzle 
Because of departures from isentropic behavior, however, stagnation pressure losses will occur These 
losses may be expressed in terms of an expansion efficiency defined as 


^(e)=7 ^ )aC ‘ Ual = ~ — — ^ 1 

* ' isentropic hy 


(39) 


where h y is the static enthalpy at the end of expansion necessary to achieve the same static pressure drop in 
an lsentropK process. In reality, the process depends upon flight conditions, but this approach is sufficient 
within the assumptions of our analysis. Equation (39) may also be written in the form: 


^ls(e) ~ 


= l -fe/r 4 ) 


h 4-hy \-(T y /T 4 ) 


(40) 


which may be solved for T IT \ 


y l 4 


jy _ Ve ~ 1 
T 4 %(e) 


+ 1 


(41) 


where we have introduced the static temperature ratio for the 
the idealized isentropic expansion process is 


expansion as \f/ e — T 5 /T 4 . Gibbs equation for 


dp( r W 

P \r-\)T (42) 

and integration of eq. (42) through the expansion gives 


P5 = Py ( ^yV-1 
P4 P4 U 4 J 


(43) 
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Substitution of eq. (41) into eq. (43) eliminates T y /T 4 


P 5 = 
P 4 


3^+1 

^s(e) 


7-1 


( 44 ) 


Using the Mach number relation between stagnation and static pressures, we can write the 
stagnation pressure ratio for the expansion in the form. 


n e = 


P0,5 _ P5_ 
POA P* 




\+A±mI 


7-1 


(45) 


Using the Mach number relation between stagnation and static temperatures, we can write the 
stagnation temperature ratio for the expansion in the form. 



To a t A\ + 


1 

y-i M i 
m 4 


But the total temperature is conserved, and eq. (46) can be written in terms of %■ 


(46) 



i+LImI 

i + ^V 


( 47 ) 


Combining eq. (47) and eq. (45), we therefore deduce 



Upon substituting for p 5 /p 4 using eq. (43), we obtain the desired final result for the stagnation 
pressure ratio: 


„ _ P0,5 _ 
7l e — — 


V'e-l , 1 


7-1 


POA 


Vs(e)Ye Ye 


(49) 


2.5.5 MHD Accelerator 

Consider, now, the MHD accelerator in which the stagnation pressure must increase. We begin with 
the Gibbs equation for an ideal (isentropic) process as applied to stagnation properties and integrate through 
the accelerator to obtain the stagnation pressure ratio in terms of the total enthalpy ratio: 


n 


a 


P 0,6 

P0, 5 


^06 

^*0,5 



(50) 


Utilizing the definition of given by eq. (6) and the definition of 7] s(a) given by eq. (8), we can 
form the following relationship between accelerator performance parameters: 


PN(a)Ps(a) — 


(*0,6“ *0,5 )* entropiA^- fys), 


actual 


(%6 “ ^,5) 


actual 


*0,5 


/ fy ),6 ~ fy ,5 ^ 


%5 


) 


isentropic 




-1 


isentropic 


(51) 


Upon substituting eq. (5 1 ) into eq. (50), we obtain an expression for the stagnation pressure ratio in 
the MHD accelerator in terms of the thermodynamic parameters of the device 


K 


a 


p 0,6 

PQ,5 


Y 

[ l + T lN^a) T ls(a))y- 1 


(52) 
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If we eliminate % (a) using eq. (20) and eliminate T 0 l using eq. (10), we obtain the desired final 
expression for the stagnation pressure ratio in the MHD accelerator 


*a = 


_ P0,6 _ 
P 0,5 


1 + fJs(a) 


1 + /. 


0 “ XVlNig) T 0,a 


l-XVNig) t oa] 


7-1 


where we have also used 7’o,5 _ ^o,4 - 


(53) 


2.5.6 Nozzle 

The irreversible losses in the nozzle are accounted for in terms of a nozzle expansion efficiency in 
the same form as defined for the cooling expansion in eq. (39): 


(factual -!%zho<\ 

nsin) ~ misentropic 


(54) 


Thus, the relationships derived for the cooling expansion directly apply to the nozzle with appropnate 
modifications to the engine station references. For instance, when the form of eq. (48) is applied to the 
nozzle expansion process, we obtain 


*n = 


P0,10 

Po, 6 


P 10 
P 6 


{Vn 


7-1 


( 55 ) 


where * s the static temperature ratio of the nozzle. In this case, however, it is advantageous to 

express the stagnation pressure ratio in terms of the nozzle static pressure ratio rather than the nozzle static 
temperature ratio (!//„). This can be done by applying the form of eq. (44) to the nozzle, solving for •//„, and 
using this result to eliminate \ff n in eq. (55). This yields the following form for K n . 


n n = 


J 7 S ( n )+( 1 


\(r-')/r 


-y 

y-i 


(56) 
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2.6 System Thrust Characteristics 


There is now sufficient information to compute the total engine thrust including the effect of 
aerodynamic losses. Although the n parameters are not truly constant over large Mach number variations 
and y can vary significantly through the engine flowpath, the following development demonstrates the 
essential performance characteristics of hypersonic airbreathing engines. 

First, we utilize the Mach number relation between stagnation and static pressures as defined by 
eq. (5) and evaluate these expressions at the inlet and exit of the engine. We then form the ratio between 
these expressions and solve for the exit Mach number M 2 '- 


Mg = 




( 

(y-O/y 

2 

fi 

P0,e P a ’ 

-l 

7-1 

V 2 ) 

^PQ,a Pe ) 



(57) 


but 


P 0,<? _ Po,l P0,2 PO, 3 POA P0,5 A), 6 P0,10 
PO,a P0,a POA P0,2 P0,3 POA P0,5 P0,6 


— 7Cp7t ( j7Cg7lf ) 7Cg7l a 7[ n 


(58) 


where n p , K d , K h , n e , n a , and n n are the n parameters for the pre-ionizer, diffuser, generator, burner, 
cooling expansion, accelerator, and nozzle, respectively. Thus, in terms of the component stagnation 
pressure ratios 


M 2 e = 


7-1 


i +*-± j #2 


\(r-i)/r 


ft p^d^g^b^e^a^n 


Pa 

Pe) 


-1 


(59) 


We now define a global stagnation pressure ratio parameter n as 


n = i 1+2— M 2 


ft p^d^g^b^e^a^n 


Pa 

Pe 


\{r-i)/r 


such that 


(60) 
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Assuming that heat transfer from the engine is negligible, the exhaust velocity is given by 
u e = M e ^fyRT e or, in terms of the exhaust stagnation temperature 


Ug — Mg, 


7 ^ 0 , 6 / 




(62) 


where we have made use of the fact that ^o,io _7 o,6- The fuel-air ratio necessary to produce the desired 
value for T Q A is given by eq. (16). 

The thrust for a hypersonic airbreathing engine is defined by 


F = (m a + m e )u e - m a u a + (p e - p a )A e 


(63) 


and the thrust per unit mass flow rate of air becomes 


F 1 

— = [(1 + f)“e ~ u a] + —{Pe ~ Pa) A , 
m n m n 


Substituting eq. (62) into eq. (64) yields the desired form: 


(64) 


F_ 

m a 


h . |2y/?7b, 6 (n-i) i 

(1 + /) V (y-l)ri M a<y RT a 


Pe A e 




Pa 
Pe J 


(65) 


where T 0 6 is obtained from eq. (22) using To ,4 ('- e -’ To,ii m ) as a parameter. The fuel specific impulse may 
also be determined as 


F _ F/w a 

sp w f f ( 66 > 

where w a and are the weight flow rates of air and fuel, respectively. Following conventional practice, 
the TSFC is defined as the ratio of fuel mass flow rate to engine thrust: 


me f 

TSFC = -*- = — — 
F F/< 


(67) 
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2.7 Ramjet Mode Mach Number Constraint 

The limitation on burner entry temperature leads directly to a restriction on the burner entry Mach 
number My That is, the flow can be decelerated only a limited amount before the static temperature at the 
burner inlet becomes too great for effective combustion, and at some particular flight Mach number, it 
becomes necessary to transition to a supersonic combustion mode (M 3 > 1). For the MHD-bypass engine 
configuration, however, the burner entry temperature is affected by the enthalpy extraction process in the 
MHD generator. Indeed, the bypassing of flow enthalpy around the burner enables the vehicle to obtain a 
higher flight Mach number while maintaining a subsonic combustion mode. 

We can easily develop a relation for the burner entry Mach number as a function of the flight Mach 
number, the enthalpy extraction ratio, and the static temperature limit. We begin with the Mach number 
relation between stagnation temperature and static temperature at the burner entrance 


^0,3 - T 3 


1 + 


( 68 ) 


Then we eliminate T 0 3 using eq. (13) while enforcing a burner entry static temperature limit (T 3 | im ) 
and find that the burner entry Mach number is given by 


M 3 = 



Tg 

T \\\m 


1 ~ OBmI 

l ~XVN(g) 


i + 


7~1 

2 




This relation implies definite operating restrictions. For example, when 


(69) 


M a < 


7-1 


^3,lim 1 X^lN(g) j 


T 

l a 


1- 




J 


(70) 


no solution is physically possible since T 3 lim would exceed the stagnation temperature of the freestream 
flow. Furthermore, when 


M a > 


7-1 


y + 1 


^ ^3,lim 1 > 

J T a 1- 


WN(g) 

Vn( 8 ) 


-1 


the flow entering the burner must remain supersonic or T 3 lim will be exceeded. 


(71) 
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3. REPRESENTATIVE PERFORMANCE CALCULATIONS 


3.1 Specific Thrust 

It is instructive to compare the performance of an MHD-bypass hypersonic airbreathing engine 
with that of a conventional ramjet system. For this purpose, Hill and Peterson’s well-known analysis for a 
nonideal ramjet represents a suitable baseline case.' 3 Indeed, their analytical results are directly recoverable 
from our development by simply eliminating all MHD interaction (i.e., setting t] N{ ) %(o) >■ 

Following Hill and Peterson, we assumed a freestream temperature of T = 220 K and used the thermo- 
dynamic properties of air throughout the engine flow path. The fuel heat of combustion was taken to be 
/=45xl0 3 kJ/kg with 100 percent combustion efficiency and a stagnation temperature limit was enforced 
at the burner exit. As in Hill and Peterson’s case, we considered constant stagnation pressure ratios of n d 
=0.7 7Ti= 0.95 and n =0.98 for the diffuser, burner, and nozzle, respectively. As previously noted, the 
calculations are limited because y is not constant throughout the engine flowpath and the stagnation 
pressure ratios are not independent of flight Mach number. These coefficients were taken as constants 
solely for the purpose of maintaining simplicity and clarity. We have optimistically assumed that 
„ „ o.9 and computed n (loss) and n a (gain) using eq. (25) and eq. (28), respectively. The 

power needed to drive the pre-ionizer is difficult to estimate at this juncture and we utilized a value of 
^=0.05 for demonstration purposes. 

The computed specific thrust is shown as a function of flight Mach number in fig. 2 forenthtilp y 
extraction ratios of % ( „)= 0, 0.25, and 0.5 and a maximum burner stagnation temperature of 3000 . e 

case without any MHD interaction corresponds directly with Hill and Peterson s results and is demons ra- 
tive of a conventional hypersonic airbreathing engine. The specific thrust in this instance peaks at t a lg t 
Mach number between 2 and 3 after which it steadily decreases to zero near M fl -8. The results with MH 
interaction demonstrate a reduction in peak performance and a shifting of the specific thrust curve to hig er 
flight Mach numbers. The loss in performance using MHD-bypass is due to increased non-, sentropic losses 
associated with the MHD energy conversion devices. However, the ability to bypass flow enthalpy aroun 
the burner permits optimization of the combustion process at higher flight Mach numbers and enables the 
production of thrust in a flight regime not obtainable with conventional engines. This capability is the 
primary advantage of the MHD-bypass concept. 
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Figure 2. Specific thrust characteristics as a function of flight Mach number for t) N , .=0, 0.25, 

and 0.5 assuming 77 j(5) =77 i(a) =0.9 and r 0 |im =3000 K. * 
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3.2 Ramjet Mode Envelope 

The combustion optimization attributes can be better appreciated, perhaps, through an inspection of 
the maximum flight Mach number which can be achieved with MHD-bypass of flow enthalpy while main- 
taining subsonic combustion conditions. This effect, shown in fig. 3, is demonstrated through solution of 
eq (71) assuming a maximum static burner entry temperature of T 3 lim =1600 K.This envelope indicates that 
the enthalpy extraction ratio in the generator must be extremely high to maintain ramjet mode operation, 
irrespective of the process efficiencies. 



Figure 3. Maximum flight Mach number for maintaining subsonic combustion conditions as a 
function of the generator enthalpy extraction ratio. Assumes a maximum static burner 
entry temperature of T 3 ljrn = 1600 K. 
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3.3 Stream Thrust Analysis 

It has not escaped the authors' attention that the present cycle analysis is greatly limited bv the 
assumption of constant specific heat and specific heal ratio throughout the engine flow path. In practice 
these parameters change dramatically. The next step in analysis sophistication is to require that these 

parameters remain constant only within each engine component. This is best done through introduction of 
the stream thrust function (Sa) defined as 


Sa = w| 1+ 


(72) 


, ivnuiTu 6 eam hrUSt approach ’ P 10neere d by Curran and Craig 14 , can be readily generalized for analysis 
o HD-bypass engines. For example, it can be shown that the specific thrust of the MHD-bypass engine 
illustrated in fig. 1 can be expressed in the form 12 (p. 175 of the reference): 


• -[( 1 + /) 5 a io -^ fl ]- — 

"In U 


R„T a 


ho 


-1 


(73) 


where 


A a Pi 0 T a «, o (74) 

By tracking the evolution of the stream thrust function through the engine, we completed preliminary 
va idation calculations for this analysis methodology. This technique can be utilized to conduct more compre- 
hensive parametric studies and to obtain more realistic estimates of MHD-bypass system performance. The 

results based on the stream thrust methodology, however, do not alter the basic scientific feasibility findings 
as deduced from the simplified model. 6 
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4. MAJOR TECHNICAL ISSUES 


The introduction of plasma/MHD technologies into the propulsion system introduces a host of complex 
technical and systems integration issues. Most prominent among these is the ionization characteristics an 
electrical characteristics of the inlet air under various flight conditions and MHD-device performance. 


4.1 Air Ionization Considerations 

The two major parameters of interest with respect to ionization quality include the electrical 
conductivity (a) and the Hall parameter (j8). In principle, the degree of MHD interaction depends upon the 
magnitude of electrical conductivity. However, it is the magnitude of 0 which can determine the optim 

design configuration. 

4.1.1 MHD Interaction Parameter 

The electrical conductivity is a critical parameter in that it determines the minimum magnetic field 
strength required for achieving meaningful levels of MHD interaction. For example lower conductivity 
implies the^need for larger, heavier, and more power-hungry electromagnets. Generally, the manifestation 
of significant MHD effects requires an interaction parameter on the order of unity. The interaction param- 
eter (0 is defined as the ratio of MHD forces to the inertial forces in the flow: 


Q = 


oB^L 

pu 


(75) 


where cris the electrical conductivity, B is the magnetic induction, L is a characteristic length, p is the flow 
mass density, and u is the flow velocity. 

Practically speaking, the product cB 2 should be as high as possible for a given geometry, as defined 
by L and for a given flight condition, as defined by the mass flux pu. Therefore, it is advantageous to 
utilize the highest achievable electrical conductivity to minimize the required magnet strength. 

In general, (Tvaries in proportion to T a lp m where a is on the order often. This implies that plasma/ 
MHD concepts work best at high flight speeds (high temperature) and high altitudes (low 
conditions are ideal for a large portion of the hypersonic flight regime. A question immediately , 
however, as to the minimum flight speed necessary for achieving a naturally occumng plasma capable 
sustaining meaningful MHD interaction at reasonable magnetic field strengths. 

To address this question, it is of interest to map out the naturally occurring interaction parameter as a 
function of flight Mach number and altitude, and calculations were performed for equihbnum air behind a 
normal shock wave. These calculations were carried out with a modified version of the NASA SP .73 


24 



chemical equilibrium code in which the plasma electrical transport properties are evaluated according to 
the methodology ot Frost. 5 > 16 Results based on the 1976 Standard Atmosphere are shown in fig 4 
assuming a magnetic field strength of 1 T and a characteristic length of 1 m 



Figure 4. Naturally occurring interaction parameter behind a normal shock as a function of 

flight Mach number and altitude. Based on 1 976 standard atmosphere and equilibrium 
air with molecular dissociation. 


Additional sets ot equilibrium air calculations have been performed both with and without cesium 
seed material. The resulting interaction parameter is shown in fig. 5 as a function of the flight Mach 
number for fixed altitudes of 1 00,000 and 1 25,000 ft. Note that for a concentration of 0.05 percent cesium 
(molar) significant interaction is available at Mach numbers <10. The dips in the curves are associated 
with molecular dissociation effects. The fundamental conclusion of these calculations is that equilibrium 

° Hnlet air 1S only P ractical at mid-to-high hypersonic Mach numbers, say 1 2 to 1 8. Application 

o MHD devices at lower flight speeds will require a power source for inducing nonequilibrium ionization 
or the air. 
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Figure 5. Interaction parameter behind a normal shock as a function of flight Mach number 
for both unseeded and cesium-seeded equilibrium air. Calculations were performed 
for fixed altitudes of 100,000 and 125,000 ft. 
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4.1.2 Hall Parameter 


In general, two distinct electric fields are induced when a moving conducting fluid interacts with an 
applied magnetic induction. The first of these is the Faraday field associated with the bulk motion of 
charged particles in an applied magnetic field and is defined by the vector £p=uxB. The second of these is 
the Hall field which is associated with the electromotive force (emf) generated by the drift velocity w of the 
charged particles and is defined by the vector %=wxB. Both of these induced fields, in addition to being 
orthogonal to each other, are orthogonal to the applied magnetic field. Furthermore, these induced fields 
give rise to two separate components of current density in the conducting fluid which are defined by the 
following component equations: 


and 


Jf - 


1 + /T 


|uxB| 


(76) 


JH = 


i + /r 


|/3|uxB| 


(77) 


Thus, the relative magnitude of the Faraday and Hall current components depends upon the 
absolute magnitude of the Hall parameter. When p exceeds unity, the Hall current becomes greater than the 
Faraday current. Thus, the magnitude of (3 determines the appropriate design configuration for any 
particular application. 

Because the Hall parameter increases dramatically with decreasing pressure (i.e., increasing 
altitude), it is important to assess the anticipated variation in /J at the flight conditions of interest. To obtain 
quantitative predictions for the the Hall parameter under high-altitude hypersonic flight, equilibrium air 
calculations were performed in which P was deduced assuming an applied magnetic field of 1 T. 

The resulting values for the Hall parameter are shown in fig. 6 as a function of the flight Mach 
number for fixed altitudes of 100,000 and 125,000 ft. One is immediately struck by the large magnitude of 
P under these conditions, and it is obvious that the results have serious implications with respect to system 
design configuration. 

For a freestream M of 10, for instance, P is ~2 at 100,000 ft and increases to ~1 at 125,000 ft. This 
means that the Hall current can be as much as 7 times larger than the Faraday current within this flight 
envelope. For an applied field of 2 T, the Hall parameter would be twice as large. Based on these results, we 
conclude that the design configuration of plasma/MHD devices must be optimized to take full advantage of 
the intrinsically large P values. 
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Freestream Mach Number 


Figure 6. Hall parameter behind a normal shock as a function of flight Mach number for 
equilibrium air. Calculations were performed for fixed altitudes of 100,000 and 
125,000 ft with B= IT. 


4.2 MHD Technology Considerations 

The technology for MHD generators has matured significantly in recent years and efficiencies are 
approaching the level needed for practical application. The technology for MHD accelerators, however, is 
comparatively less developed and there are lingering uncertainties with respect to achievable performance. 

A small number of experimental efforts have been conducted over the past 40 yr, but the results 
have not been conclusively positive.* 7 - 24 In many cases, the resulting velocities were much lower than 
theoretically anticipated. In fact, it has been noted that the observed flow accelerations m these experiments 
could be associated purely with a thermal heating mechanism rather than true MHD interaction. 

In a supersonic flow, for instance, the temperature in the boundary layer near the insulated sidewalls 
of an MHD device can exceed the temperature in the inviscid core, and the electrical conductivity in the 
boundary layer will be correspondingly higher. Furthermore, the velocity in the boundaiy layer is ower 
than in the inviscid core, and the back emf will be lower. One is therefore led to conclude that there will be 
less resistance to electric current flow in the boundary layer and that severe leakage currents will occur. 

This issue as well as other practical concerns such as maximum axial voltage gradient and 
generator-accelerator coupling imply that serious research and development must be earned out on the 
technology of MHD accelerators before the practical feasibility of an MHD-bypass propulsion system can 

be assessed. 
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It should also be noted that the development of MHD systems for hypersonic aircraft is currently 
limited by available enabling technologies for flight-weight magnets and power conditioning units. The 
primary limitations are related to available magnet field strengths and materials, as well as the intrinsically 
high weights of magnet materials. In general, superconducting magnet systems are required so that power 
consumption does not become impractical; unfortunately, this entails the utilization of liquid helium 
cooling. As an alternative, it has been suggested that it might be possible to employ high-temperature 
superconductors (HTSC) which can utilize liquid hydrogen as the coolant. This would be desirable for 
hydrogen-fueled aircraft, however, fabrication techniques and critical current densities for HTSC materials 
are very limited in their current state of development. Clearly, further research and development are needed 
on flight-weight magnet enabling technologies for MHD-based aerospace systems. 
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5. CONCLUSIONS 


It has not escaped the authors’ attention that our simplified thermodynamic cycle analysis is greatly 
limited by the assumption of constant stagnation pressure ratios for the various engine components, 
nevertheless, the analysis does reveal the essential operational benefits of MHD-bypass engines and 
demonstrates fundamental scientific feasibility. Clearly, bypassing energy around the burner extends the 
operational Mach number envelope of conventional airbreathing engines, but it is important to note that 
system performance is extremely sensitive to non-isentropic losses in the MHD devices and that any 
favorable operational characteristics will disappear if these losses become too large. We conclude that 
MHD-bypass systems hold significant promise for extending the effective operating range of conventional 
hypersonic airbreathing engines and believe that the present results justify a more in-depth analysis in 
which the stagnation pressure ratios become dependent on flight Mach number and the combustion process 
is modeled more accurately. 

In this work, calculations indicated that equilibrium ionization of inlet air is only practical at 
mid-to-high hypersonic Mach numbers, say 12 to 18, and at altitudes exceeding 100,000 ft. The 
introduction of ionization seed can relax the Mach number requirement to about 10, but at the expense of 
increased propellant mass fraction. Calculations also indicate that the Hall Parameter (p) will become 
significant at these same flight conditions, and it is clear that the design configuration of MHD devices 
must take into consideration intrinsically high P values. 

A review of MHD technology revealed that MHD generator technology is relatively mature whereas 
MHD accelerator technology suffers from lingering uncertainties with respect to achievable performance. 
The major concern appears to be the effectiveness of Lorentz force acceleration in MHD devices. In past 
experiments, for instance, it is suspected that observed flow accelerations may be largely the result of a 
thermal heating mechanism associated with current leakage in the boundary layer. Further experimental 
research is needed to resolve this critical issue. 
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